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We introduce a one-dimensional system of fermionic atoms in an optical lattice whose phase 
diagram includes topological states of different symmetry classes. These states can be identified by 
their zero-energy edge modes which are Majorana fermions. We propose several universal methods of 
detecting the Majorana edge states, based on their genuine features: zero-energy, localized character 
of the wave functions, and induced non-local fermionic correlations. 



There is at present growing interest in topological 
phases of matter - fermionic insulators and superconduc- 
tors with gapped excitation spectrum and nonlocal topo- 
logical order (TO). TO is characterized by a topological 
invariant taking discrete values [T]-[6]. A striking conse- 
quence of TO is the existence of robust zero-energy modes 
localized at defects, edges, and interfaces between differ- 
ent topological phases [7H9]. The number of these states 
and their properties are directly linked to TO and, there- 
fore, can be viewed as a probe of the topological state 
providing access to the corresponding phase diagram. A 
striking example are zero-energy Majorana fermions, as 
discussed in the context of topological quantum comput- 
ing p!QHT2] . Being widely discussed in both high energy 
and condensed matter theory, the observation of Majo- 
rana fermions has so far been elusive, and thus remains 
an outstanding challenge for both experiment and theory. 

Most of the developments and proposals in realizations 
and studies of topological states of matter in general, and 
Majorana fermions in particular, have been related to 
condensed matter systems (see e.g. p!3HT7j ). Many-body 
systems of cold atoms and molecules, on the other hand, 
promise not only a clean realization of these phases due 
to the control available there p!8H23] , but also a plethora 
of new measurement techniques [2T , as illustrated by the 
recent achievement of single site addressing and detection 
in optical lattices [25l[2^. In light of these advances, we 
study below the signatures of Majorana fermions, and 
the related measurement of topological phases and phase 
transitions by available experimental techniques, includ- 
ing time of flight [27 and probe spectroscopy. This com- 
plements the analysis of the topological invariant in the 
bulk as proposed in j28ti3Q] . We illustrate our ideas in 
the context of two examples. We first propose an atomic 
"zig-zag Kitaev chain" in an optical lattice (c.f. Fig. 
la)): this provides a model with a rich topological phase 
diagram, allowing different topological states in two sym- 
metry classes supporting Majorana edge modes, as a test 
ground for various detection techniques. The presence of 
TO and the topological phase transitions can be studied 
by detecting the related Majorana edge states. We dis- 
cuss several universal methods to detect them via their 
main features, (i) zero-energy, (ii) localization near the 
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FIG. 1. (Color online) (a) A zig-zag chain of spinless fermonic 
atoms according to the Hamiltonian H. This geometry allows 
for an experimental control over the relative strength of NN 
(red) and NNN couplings (grey). C2j-i and C2j denote Ma- 
jorana operators on lattice sites j — 1 to L (see text), b) 
Preparation of the initial state of Majorana fermions by adi- 
abatically cutting the closed chain. 



edge and (iii) induced non-local fermionic correlations. 
Finally, we extend these ideas to the case of an atomic 
Px + Wy superfluid. 

Model of a zig-zag Kitaev chain: We consider single- 
component fermions on a finite chain (one-dimensional 
(Id) lattice) of size L with Hamiltonian H = Hi -\- H2 
H^. Here, 
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where and a J are fermionic operators {i = 1, . . . ,L), 
Ja > and Ac, = |Ac,|e'^" are nearest- (NN, a = 1) 
and next-to-nearest-neighbour (NNN, a = 2) hopping 
and pairing amplitudes, respectively, and fi is the chem- 
ical potential. H can be viewed as a Hamiltonian of 
two coupled Kitaev chains [31 (with odd or even sites). 
In Fig. [1^) we represent the chain as a zig-zag, which 
leads naturally to a cold atom implementation with opti- 
cal lattices [22 as the relative strength of NN and NNN 
amplitudes can be controlled by changing the zig-zag ge- 
ometry. The pairing terms in are obtained by a Ra- 
man induced dissociation of Cooper pairs (or Feshbach 
molecules) forming an atomic BCS reservoir, with /i a 
Raman detuning. 

The symmetries of H and the corresponding sym- 
metry class [1 depend crucially on the relative phase 
(j) = — (j)2 of the pairing amplitudes (one of the phases, 
say, can be gauged away by redefining the operators 
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FIG. 2. (Color online) a) Phase diagram corresponding to 
the Hamiltonian if for /x = 0, |Ac| = J^e'^. b) The two 
lowest eigenvalues £1,2 of if as a function of the ratio J2/J1: 
If = 0, TT (mod 27r) (BDI class) there are one or two such 
modes (b), for all other values of cj) (D class), there is at most 
one zero mode (c). 



aj). For = 0, TT (mod27r) it is in the class BDI with 
time-reversal, particle-hole, and chiral symmetries, for all 
other values of (j) it is class D with only particle-hole sym- 
metry [38 . Both classes allow topologically nontrivial 
states characterized by the integer-valued winding num- 
ber G Z in the class BDI and by the Chern parity 
number Pc = 0, 1 G Z2 in the class D. Since we dis- 
cuss the topological phase diagram of H in full generality 
elsewhere, we summarize here only the relevant features 
for our detection schemes. The phase diagram contains 
states with u = 0,±1,±2 for the BDI class, and with 
Pc = 0,1 for the D class. Fig.|2|i) summarizes results for 
the simplest case ja = and = Ja^ where the phase 
diagram is characterized by the relative phase (j) and the 
ratio J2/J1' We find that u = Pc = I ioi Ji > else 
one has Pc = and u = -\-2 or —2 for ^ = and tt, 
respectively. Therefore, for both symmetry classes, the 
point Ji = J2 corresponds to a topological phase transi- 
tion, where the bulk excitations become gapless [38 . 

The presence of a nontrivial TO in the bulk leads to 
zero-energy modes at the interfaces between states with 
different TO. The number of these states is related to 
the difference of the corresponding topological invariants 
[HI [9], which in our case are Majorana fermions. Using 
the hermitian Majorana representation, where C2j-i = 
Oj + aj and C2j = {—i){aj — aj) obey {cj^ci} = 25 ji, 
the quadratic Hamiltonian reads H = i ^ AjiCjCi with 
real A = —A^. Diagonalizing A for a finite chain we 
identify the zero-energy modes and associated Majorana 
edge- mode operator = "^k^Ek^k with (real) coeffi- 
cients VEk localized at the interface with some localiza- 
tion length lioc [To be precise, the energy of these modes 
~ exp(— i^//ioc) approaches zero for L 00.] 

For the BDI class we have in total two such modes 
for Ji > J2: 7l = ci and 7^^ = C2l for the left and 
the right edge, respectively, and four modes for Ji < 
J2. The two additional zero modes are exponentially 
decaying inside the bulk (see Figj2]3) and [38 ). For 
the D class there are two zero-energy modes for Ji > J2, 
which decay exponentially inside the chain, and no such 
modes for Ji < J2 (Fig. (2)3)). 



The presence of Majorana zero-energy modes leads to 
a degeneracy of the ground state: The dimension of the 
zero-energy subspace is 2^^^/^], where Nm is the number 
of the Majorana modes (which is always even in our case) 
and [x] is the integer part of x. A state in this subspace 
generically has correlations between Majorana operators 
belonging to different edges of the system, resulting in 
nonlocal fermionic correlations, in contrast to local corre- 
lations in the bulk. For the case with only two Majorana 
modes {u = 1 or Pc = 1), these nonlocal correlations are 
related to the fermionic parity: {G± \ ^lIr \G±) = ±i for 
the ground state \G±) with even (+) and odd (— ) number 
of fermions, in full analogy to Ref. [3T. The number of 
Majorana zero-energy edge modes is determined by the 
topological invariant in the bulk, and by detecting them 
we can access the full topological phase diagram. 

Preparation. Before turning to the detection of the Ma- 
jorana modes, we briefly comment on preparation proce- 
dures 1^ 1^ ^5 . The preparation of a unique state is 
complicated by the fact that these modes are associated 
with a zero energy subspace that is insensitive to local 
perturbations. Thus, starting from a thermal state that 
has generically no long-range correlations, the available 
local processes cannot create the latter, and we end up in 
a completely mixed state in the Majorana subspace. This 
problem can be overcome by starting with a closed quan- 
tum chain in a fully paired (and, hence, unique) state 
corresponding to the ground state of a Hamiltonian (see 
Fig. lb)). This state has an even parity (all particles 
are paired) and only local correlations between the sites. 
Making use of local addressing [25l[26] we can "cut" the 
chain: We adiabatically empty a part of the chain (with 
d sites) by locally raising the chemical potential [11 . As 
a result, we obtain an open Kitaev chain in the even 
parity ground state |G+) with a non-local correlation 
(G+l jLlR = i across the chain - a unique feature 
of the presence of Majorana fermions. For the zig-zag 
we start with two decoupled subchains (with even and 
odd number of sites) corresponding to Ai = Ji = and 
perform the "cut". This leads to two decoupled chains 
with even parity supporting four Majorana modes. This 
condition is preserved if we turn on the couplings Ji and 
Ai between the subchains until we reach the point of the 
phase transition. Then, we end up in a state of even par- 
ity for the whole chain supporting two Majorana modes 
[38 . Note that the cutting procedure avoids the need 
to establish the edge-edge correlations via a sequential 
process, which takes polynomial time due to the Lieb- 
Robinson bound (35] [36] . 

Time- of- Flight (TOP) imaging. Non-local correla- 
tions can be detected in TOF imaging, as illustrated 
in Fig. ^i). We consider atoms of mass m released at 
time t = from lattice sites Rj = ja with lattice spac- 
ing a. At time t mo? /h (far fleld approximation) 
the atomic density distribution at the detector is given 
by {n{x)) - Y.j,j' e^'<^^-^o')^^^^''\a]aj>) revealing the 



3 



a) Ri = aj d)2 



4 J2/J1 1 



.b) 


■-bulk 

— Majoranas 


/ 


c) 


■-H = 0.5 
— |i= 1.5 










\ 

\A 


■■■^ = 2.5 


a/ 


AA A A 

N \/ \ A A / \A' 




A V\ /V A 


-80 


X 


80 


-80 


X 


80 



FIG. 3. (Color online) a) Principle of TOF. The signal can 
be enhanced by excluding the bulk (grey shade), b) TOF for 
the ideal Kitaev chain: The presence of the Majorana modes 
(solid) leads to oscillations of the signal S{x) (see main text) 
which are absent in the bulk (dashed), and vanish if we ap- 
proach the transition point at /x = 2, see c). d) The amplitude 
of S(x) allows to determine the number of Majorana modes 
(see main text). 
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FIG. 4. (Color online) a) Spectroscopic detection of Majorana 
modes in a two-band lattice (see text), b) The momentum 
distribution of atoms in the upper band reveals the presence 
(left) or absence (right) of Majorana modes (see text), c) 
From the strength of the signal {uk) at resonance (here: k = 
7r/2) we can deduce the number of Majorana modes. 



initial correlations of atoms in the lattice. Therefore, the 
long-range Majorana correlations Rj — Rj' ~ Lm (dis- 
tance between the edge modes) will result in rapid oscil- 
lations of {n{x)) as compared to slow oscillations origi- 
nating from short-range bulk correlations {Rj — Rj' ^ a). 
The contrast of the Majorana signal can be enhanced by 
local addressing: we shade all but A^^ sites adjacent to 
each of the two edges (A^5 = 2 in Fig. [3|i)), thus measur- 
ing {n{x)) = 2Ni)S{x)^ where we have normalized by the 
number of sites contributing to the signal, 2Ni). 

For Ja = AcK, J2 = /i = (ideal Kitaev chain), we find 

(n(x))+ Mhulk{x) + m{d)J\fedge{x), 

where Mi,uik{x) = 2msm\^) + cos(^) stems 
from the atoms in the bulk only and the two edge 
sites separated by Lm = ctd give rise to Medge{x) = 
cos(^^^^(i)/2 multiplied by the Majorana correlation 
m{d) — —i^^LlR] [38]. Therefore, the second rapidly 
oscillating term provides information on the presence of 
the long-range Majorana correlations related to the oc- 
cupation/purity of the Majorana subspace. If the initial 
state of the chain is prepared as described above, then 
m{d) = — z (G+l 71,7^^ = 1, and the amplitude of 

this terms reaches its maximum. The results of the cal- 
culations for the ideal chain are shown in Fig. jsj)): The 
presence of the Majorana modes leads to oscillations that 
are absent in the bulk. By ramping the chemical poten- 
tial, we can detect the location of the phase transition, 
since the oscillations disappear when crossing the tran- 
sition point at /i = 2 (see Fig. [sj:)). Further, from the 
amplitude A of the signal S{x) we can deduce the number 
of Majorana modes, as depicted in Fig. [sji): The ampli- 
tude for four zero modes (J2 > J\) drops by a factor of 
two upon reaching the transition point J2 = Ji- 

Spectroscopy. The combination of a spectroscopic 



setup with TOF allows to measure not only the energy of 
the Majorana states but also their wave functions. Con- 
sider a Id lattice with two bands (modes aj {bj) in the up- 
per (lower) band) separated by an energy difference A£^5. 
In the lower band we realize the Hamiltonian while 
the fermions can hop freely with the hopping amplitude j 
in the upper band. The full Hamiltonian reads Hq = H-\- 
Hup, where Hup = j Yd^^{h\hij^i^h.c.)^ /\Ei, Yd=i ^l^i- 
We then couple upper and lower bands on the first and 
the last sites of an open chain via a time-dependent per- 
turbation V{t) = Vo{b\ai^b{aL)e-'^^ -\-h.c. (see Fig.ji^) 
and measure the momentum distribution in the upper 
band as a function of the frequency ^7, e.g., via TOF 
imaging. Assuming an ideal Kitaev chain in the lower 
band, the upper band momentum distribution is [38] 

K(t,(^)) = Fo'M-i(no + ni), (1) 
^0 = Aoft{sk - ^^), ni = ft{sk + 2J - Qh), 

where we have introduced the notation Aq = {m-^{d) — 
m-{d)) sin^ | (L - 1) + A_ and (x) = 4 sin^ [xt /2h) /x'^ . 
Further, = AEi^ + 2j cos k is the dispersion in the up- 
per band and m±{d) denote the initial occupation of the 
even/odd parity subspace in the lower band. The contri- 
bution no corresponds to a parity flip in the lower band 
(no energy cost) and the creation of an excitation in the 
upper band (energy Sk)- Thus, the absorption peak at 
hft = Sk in no indicates the ground state degeneracy. 
The edge-edge Majorana correlation length is encoded 
in the oscillation period. In contrast, the term rii re- 
sults from creation of excitations in both upper (energy 
Ek) and lower (energy 2J) band. The corresponding ab- 
sorption peak is located at hVt = 2J ^ providing the 
direct measurement of the pairing energy gap (2 J here) 
from the distance of two peaks. Thus, the topologically 
non-trivial phase leads to a clear signal at Vth — 2J -\- ja 
(for /i 7^ one has to replace HQ hVt — /i, cf. Fig.|4]3)). 
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FIG. 5. (Color online) Measurement of the localization length 
of the Majorana mode via spectroscopy (see main text). 



Note that this setup allows to detect the presence of the 
zero energy Majorana subspace independently of its pu- 
rity. Further, the number of Majorana modes is encoded 
in the strength of the resonance signal for a fixed mo- 
mentum k. In Fig. ^) we present the results for the 
Hamiltonian H with = |Aq,|, ^ = 0, where the pertur- 
bation is on sites 1, 2, L — 1 and L. Starting from a state 
with four Majorana modes, we reduce the ratio J2/J1. 
At the transition point, the amplitude of the resonance 
signal (n/c=7^/2) drops to half its magnitude indicating the 
disappearance of two of the Majorana modes. 

With a slight modification of the above setup one can 
also infer information about the localization length of the 
Majorana zero modes. Namely, we apply the perturba- 
tion to the first ul sites: V^^^)(t) = Vq YTjLi a]hje'^^ + 
h.c. The corresponding momentum distribution in the 
upper band is now 



L-l 



E 



(2) 



where Vji are the coefficients of the Bogoliubov trans- 
formation to the quasiparticle operators clj, {v labels the 
quasiparticle modes with energies Ey) in the lower band, 
ai = ^jiuiiyCLiy — v*^al), which diagonalizes H. Let us 
now consider k = and the external frequency being 
in resonance with the lowest energy {u = Majorana) 
level, Q = Qr = A£^5 + ^k=0' In this case, only the 

Majorana mode is excited in the lower band, and one 

2 

has (n/c=o(^r, ^)) ~ Yl^Li^jO = M{nL). The local- 
ized character of the Majorana mode can now be estab- 
lished by looking at the dependence of (n/e=o(^r, ^)) on 
the number of sites ul affected by the perturbation. To 
be more specific, the absolute change of M.{nL) when ul 
is increased by one, AA^(nL) = \M.{nL + 1) — M.ij^L)\^ 
is expected to be AA^(nL > 1) = for the ideal Ki- 
taev chain and AA^(nL) ~ exp(— n^/Zioc) for a generic 
one. This is reflected in a comparison of the normalized 
quantity AA^(ni,)/AA^(l) (red circles) with the normal- 
ized wave function of the Majorana mode, |vjo|^ / l^2o|^ 
(dashed) as shown in Fig. [s] 

Probing in 2d. To illustrate the universality of the pro- 
posed methods, we apply them to a 2d Px-\-ipy superfluid 
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FIG. 6. (Color online) a) Spatial distribution of the Majorana 
modes in the vicinity of two vortices imprinted on a + ipy- 
superfluid for /x = 2 J on a 32 x 32 lattice, b) A shade (grey) 
is introduced to enhance the signal contrast for the TOF (c) 
and spectroscopic (d) detections {AEb = J,j — O.IJ). 



[321 [37] with Hamiltonian Hs = H^opp -r -l-l pair ^ ^^/i, 
with a NN hopping H^oppi a pairing term Hpair = 

5^x,y ^x,ytti«J. + h.c. where Ax,y = Ax(^y,x+ei + 

^(^y^x+es)^ and Hjj^ = — /iXlx^i^x- The fermionic op- 
erators ttx are defined on a 2d (square) lattice and the 
vectors ei^2 are two independent elementary lattice trans- 
lations. We introduce two vortices (each with vorticity 1) 
via a position dependent order parameter Ax = \A\e'^^^^^ , 
which can be written onto the superfluid by phase im- 
printing in the Raman process [24 . 

For <C 4J we have a topological ground state sup- 
porting two Majorana modes located in the vicinity of the 
two vortices (cf. Fig. |6]a)). The comparison of the TOF 
signal with and without vortices (Fig.joj^)) shows that the 
presence of the Majorana modes leads to an oscillatory 
behavior of the density {n{x)) at the detector. Again, 
the contrast of the signal is enhanced by shading all but 
a small region near the vortices (Fig.[6]3)). To detect the 
Majorana modes in the spectroscopic setup, we apply an 



external perturbation V{t) = Vq ^ 



+ h.c, 



where A includes the lattice sites inside a w x w small 
regions around the vortices (here: w = 3). The result- 
ing momentum distribution in the upper band {uk^^ky) 
for kx^y = ^ is depicted in Fig. l6]i). The presence of 
the Majorana modes leads to a clear signal at fth/J = 
£JL jl/J ^ 1.19 that is absent in a vortex-free setting. 

32 '32 

In summary, we have analyzed various scenarios cre- 
ating topological states of matter with cold atoms, and 
their reliable detection via time of flight and probe spec- 
troscopy of zero energy Majorana edge modes. 
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SUPPLEMENTARY MATERIAL 
Symmetries and phase diagram of the model 

Following the general tenfold classification scheme 
[U |3i , the topological class to which a given Hamiltonian 
belongs to, is determined by the invariance properties 
of the Hamiltonian under time-reversal, particle-hole (or 
charge conjugation), and chiral (or sublattice) symme- 
try. This class specifies then, whether the states with a 
nontrivial topological order could exist and provides the 
corresponding topological invariant to distinguish them. 

For a translationally invariant Id spinless Bogoliubov- 
de Gennes lattice Hamiltonian in the quasimomentum 
basis 

keBZ 

where 

with a/e = L~^^'^^^e^^^ aj^ the invariances are equiva- 
lent to the following conditions: 

U].1-C_jJJt = T-Lk, 
for the time-reversal, 

for the particle- hole, and 

for the chiral symmetries, respectively, where Ut^ Uc^ 

and S are unitary matrices satisfying UtU^ = ±1, 

UcU^ = ±1, and S2 = 1. 

For the Hamiltonian H from the main text one has 
= —(^1 cos/c + J2 cos2/c)— /i/2 and = —i{Ai sin/c+ 

A2 sin2/c), such that the matrix Hk is 



Hk 



hk • a. 



where a = {(Jx^cry^az} is the vector of Pauli matrices 
[Im(A/e), Re(A/c), ^/c]. The excitation energies 

hk 



and h^ 



are Ek 



For generic Ai = |Ai|e*^i and A2 = |A2|e'^^ only 
the particle-hole symmetry condition is fulfilled with 
Uc = (as it should be for a general Bogoliubov-de 
Gennes Hamiltonian), and, therefore, the Hamiltonian 
H belongs to the class D. For ^ = — 02 = or tt 



(mod27r), that is Ak 



-le 



2(6 I All sin fe+ |A2|sin2fc) 



with £ = ex.p{i(j)) = ±1, the conditions for the time- 
reversal and chiral symmetries are also satisfied with 



Ut = diag(e-^^2^e^^2) ^nd S = a^Ur, and the Hamilto- 
nian H belongs to the chiral BDI class. Note that the 
phase ^2 can be gauged away by ak e^^'^^/s^^ such 
that Im(A/e) = and the vector hk belongs to the yz- 
plane for all k. Geometrically speaking, for the BDI class 
the vectors hk for all k are in the same plane, which is 
the yz-pldine for 02 = 0. 

For the chiral BDI class in Id, different topological 
states are classified by the integer valued winding number 
z/ G Z defined in the following way (see, for example, 
[7]): By gauging away the phase (j)2 and using a unitary 
transformation U = exp(z7r<j2y/4), the Hamiltonian can 
be transformed into a canonical form 



nk = 



Qk 



ml 











with Qk = ^k Ak and (fq^^ being the phase of qk. Then 



dk 



Alternatively, the winding number can be defined us- 
ing the unit vector fik = hk/ 
is gauged away) 



hk 



in the ?/z-plane (^2 

1 r 

= -- / Cx ■ {fik X dknk)dk G Z, 

as number of times the vector fik wraps around the unit 
circle when k goes around the Brioullin zone (here Cx is 
the unit vector along the x-axis). Note that as long as 

hk ^ (or, equivalently, \qk\ / 0), which corresponds 
to the gapped spectrum of the Hamiltonian, the wind- 
ing number is well-defined and, therefore, can change its 
value only when the gap closes. This gap-closing condi- 
tion is the necessary one for the topological phase tran- 
sition. 

For |Ac^| = Ja and /i = 0, the gap closing condition 
corresponds to Ji = J2, and straightforward calcula- 
tions give V 

3sign (A 



sio 



and z^0=7r = 



-_J2 



+ J2 



1 



For Ji = J2 we indeed have a 

topological phase transition between two different states. 
For a finite chain, the winding number manifests itself in 
the number of zero-energy modes localized near the chain 
edges, see Figs. [7] a) and[7|b). When one approaches 
the topological phase transition point, say, from the side 
J2 > Ji, two out of four zero-energy modes start to pro- 
liferate inside the bulk of the chain and become gapped 
on the other side of the transition; see Fig. [7|c). 

The topological phase diagram in this chiral class for 
the case |Ai| 7^ Ji, but still A2 = J2 and /i = 0, 
is shown in Fig. [s] In this case hk = (0,Aisin/c + 
J2 sin 2/c, Ji cos k -\- J2 cos 2k), where we assume ^2 = 0, 
and, hence, real Ai. The boundaries between different 
topological phases results from the condition hk =0. 
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FIG. 7. Spatial distribution of the zero energy modes of H 
for Ja = I Ac I and = 0. a) If J2 < Ji we have two zero 
energy modes that are located at the boundary, b) If J2 > 
Ji, two additional zero energy modes with an exponential 
decay appear, c) These modes extend more and more over 
the lattice when we approach the transition point J2 — Ji- 




FIG. 8. Topological phase diagram in the chiral class for the 
case I All 7^ Ji but still A2 = J2 and /x = 0. We see that 
we can obtain regions with = (white), 1 (striped) and 2 
(grey). 



For the D class [when (j) ^ 0, 7r(mod 27r)], there are 
only two different classes of the topological states char- 
acterized by the Chern parity number Pc G Z2. To 
define this topological invariant [4 , one has to extend 
the Hamiltonian 1-Lk and thus hk from a one-dimensional 
Brillouin zone k G [— tt, tt] ^ 5*^ (topologically equivalent 
to the circle 6*^) to a two-dimensional {2d) Brillouin zone 
/c, t G [— TT, Ti] xS^ = (topologically equivalent to 
a two-dimensional torus T^), 1-Lk ^k,t = hk^t ' ^5 such 
that the resulting Hamiltonian is gapped and belongs to 
the D class in 2d: cFxW_^ -t^x — —^k,t' The unit vec- 
tor fik^t = hk^t/\hk,t\ maps then the 2d Brillouin zone 
into a 2d sphere 5^. Different topological classes of such 
mappings are distinguished by the integer-valued Chern 
number 



C 



dkdt 
47r 



« 
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FIG. 9. Spatial distribution of the zero modes for = 7r/3 
and J2/J1 = 0.2. There are two zero modes (left and right) 
which have are linear combinations of only even or only odd 
Majorana operators living on a site j. 



It turns out [4 , that the parity Pc of C does not depend 
on the chosen extension of the Hamiltonian, and, there- 
fore, provides the topological invariant for Id Hamiltoni- 
ans in the D class. (We set Pc = 1 for odd C and Pc = 
for even C. Topologically nontrivial states correspond to 
Pc = 1.) 

The extension can be constructed using a hermitian 
matrix f}{k,t) = ()i(/c, t)cri, which has a gapped spec- 
trum and depends on two parameters k G [— tt, tt] and 
t G [0,7r] such that 5(A:,0) = Hk and ^(/c,7r) = IhkWz- 
Note that ^{k^t) interpolates between the initial Hamil- 
tonian Hk and the "trivial" one \hk\(Jz, both from the D 
class (for t 7^ or tt, ^(/c, t) does not necessarily belong to 
the D class in Id). Such a matrix always exists because 
it belongs to the A class of general hermitian matrices 
(with no symmetries) with a gapped spectrum that is 
topologically trivial. Then the matrix 



n 



k,t 



(l(A:,7r + t), fort G [-tt, 0] 

-aii)^{-k,7r - t)ai for t G [0,7r] 



is in the D class and provides the desired extension. The 
interpolation i){k,t) can be obtained, for example, in the 
following way: ()i(/c, t) for t G [0, 7r/2] corresponds to a ro- 
tation of hk = ()i(/c,0) to the ?/-axes, such that [)i(/c,7r/2) 
is parallel to the ?/-axes for all /c, and then ^(/c,t) for 
t G [7r/2,7r] describes the rotation of ()i(/c,7r/2) to the 
2:-axis. 

Following this strategy for the case |Aq,| = with 
/i = we find C = for J2 > Ji and C = sign(sin^) 
for Ji > J2, where (j) = ^1 — ^2- The point Ji = J2 
corresponds to the topological phase transition, at which 
the gap vanishes, see Fig. 2c). 



Preparation of Majorana fermions in Id 

We present here some details on preparation of a 
unique Majorana fermion state in an open chain. The 
idea can be illustrated on ideal Kitaev chain (A2 = J2 = 
and Ai = Ji). Starting from the ground state of 
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FIG. 10. a) Evolution of energy spectrum of the Hamiltonian 
= Hk + M(t)a^aL|j=i,L for an adiabatic ramp of /i. b) 
Evolution of the occupation and Majorana correlation func- 
tions during the adiabatic ramp. 



a closed chain, we cut it adiabatically by ramping the 
chemical potential on the site L, 

= H -\- /j.{t)a]^aL' 

This process affects only the four Majorana modes 
C2L-27 C2L-I7 ci, and the corresponding instanta- 
neous eigenvalues are and 2/i for the odd parity sec- 
tor, and /i ± + /i^ for the even one, such that there 
are two degenerate ground states of the Kitaev chain for 
/i ^ 00. Because a superposition of fermionic states with 
different parity is forbidden by superselection rules, the 
corresponding eigenstates never mix during the adiabatic 
ramp (cf. Fig. plo^)), and, hence, the final state of the 
open chain has even parity. In Fig. [IOJd) we depict the 
evolution of the Majorana correlation functions between 
nearest- {{iciC2L)^ solid) and next-to- nearest {{iciC2L-i)^ 
dashed) neighbors, as well as the occupation of the site 
L {{ul)^ dotted) during the adiabatic ramp. 

For the "zig-zag" chain we can start with two decou- 
pled closed chains (with even and odd lattice sites) when 
only A 2 and J2 are nonzero but Ai = Ji = 0. We then 
cut them as described above into two open chains each 
with two edge Majorana modes and in the even fermionic 
parity state with all fermions being paired. Now we can 
switch Ji and Ai on adiabatically to create a single quan- 
tum chain with even parity. Note that during this process 
there is no single-particle transition between the two sub- 
chains (with even and odd sites) because such a transfer 
would create single-particle excitations in both subchains 
which cannot occur in an adiabatic process due to the 
conservation of energy. Thus, each of the subchains re- 
mains in the even parity state. 

Further increase of NN amplitudes results in a topo- 
logical phase transition with closing the single-particle 
excitation gap allowing for a single-particle exchange be- 
tween the subchains. Then, only the (even) parity of the 
entire chain is preserved. 

Note that a similar idea can be used to create the initial 



Majorana state in the dissipative setting by adiabatically 
emptying one site [23^. 

TOF imaging 

The calculation of the density distribu- 
tion at the detector in a TOF experiment, 
{n{x)) = 5:^.^., e2^^(^^-^^')^/(^^)(a]a,v), where 

= 2mL/ht^ requires the knowledge of the first 
moments {o^jdj') of the state in the lattice before 
the trapping potential is switched off. For the ideal 
Kitaev chain, the non-vanishing first moments in the 
Majorana language are {c2jC2k)± = {c2j+iC2k+i)± = hj, 
{c2jC2k+i)± = -{c2j+iC2k)± = -ihj, and the edge-edge 
correlations are given by (c2lCi)± = im±{d) [m{d) = 
if we consider the bulk only] with d being the distance 
between the two Majorana edge modes. When only 
Nij sites on the right and on the left sides of the chain 
contribute to the signal (other sites are covered with a 
shade), the atom density distribution at the detector 
reads 



(n(x)) 



27V5sin^(^) 



cos(^) + ^m{d) cos(^d) 



For the general case, we obtain the correlation matrix 
{ojaj') via a numerical diagonalization of the Hamilto- 
nian H. 



Spectroscopic detection of Majorana modes 

We give the derivation of Eq. (2) for the ideal Kitaev 
chain. To this end we have to consider the effects of the 
perturbation V{t) = {b\ai + b]^aL)e~'^^^ + h.c. on the 
Hamiltonian H = Hk + A£^5^^^6j6j, where Hk = 

i2Ji ^f^i C2iC2i-\-i is the ideal Kitaev chain Hamilto- 
nian in the Majorana basis. Note that the Hamilto- 
nian Hk is diagonal in the Bogoliubov quasiparticle ba- 



sis di = {a] - ai ^ a]^^ - a^+i)/2 = (c2^+i + ic2i)/2, 
Hk = '^Ji'^f^i djdi, and its ground state is two- fold 
degenerate, with the ground states |^+) and |^-) hav- 
ing different fermionic parity. 

We can now rewrite the external perturbation in terms 
of quasiparticle operators di and Majorana operators 7^, 



usmg ai 



ai)/2, and aL = (^l,_i + 



dh-i — ^7i?)/2. Choosing the even parity ground state 
|^+) as an initial state and using 7l|^±) = 
7i?|^±) = a straightforward application of the 

standard time-dependent perturbation theory results in 
Eqs. (2) and (3) in the main text. If we go away from 
the ideal case, the ground state and the excitations can 
be determined numerically. 



